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ABSTRACT 

We present photometric observations of two transits in the WASP-50 planetary system, ob- 
tained using the ESO New Technology Telescope and the defocussed-photometry technique. 
The rms scatters for the two datasets are 258 and 211 ppm with a cadence of 170 to 200 s, 
setting a new record for ground-based photometric observations of a point source. The data 
were modelled and fitted using the PRISM and GEMC codes, and the physical properties of the 
system calculated. We find the mass and radius of the hot star to be 0.861 ± 0.057 M Q and 
0.855 ±0.019 R©, respectively. For the planet we find a mass of 1.437 ±0.068 Mj up , a radius 
of 1.138 ± 0.026 Rjup and a density of 0.911 ± 0.033 pj up - These values are consistent with 
but more precise than those found in the literature. We also obtain a new orbital ephemeris for 
the system: T a = BJD/TDB 2 455 558.61237(20) + 1.9550938(13) x E. 

Key words: planetary systems — stars: fundamental parameters — stars: individual: WASP- 
50 — techniques: photometric 



1 INTRODUCTION 

Mayor & Queloz (1995) discovered the first planet orbiting a nor- 
mal star outside our own solar system, 5 1 Peg, from radial veloc- 
ity (RV) measurements. The first transiting planet was found using 
photometric observations of a system, HD 209458, already known 
from RV measurements to host a planet (Henry et al. 2000; Char- 
bonneau et al. 2000). OGLE-TR-56 was the first planet discovered 
from its transits (Udalski et al. 2002; Konacki et al. 2003). The 
transit detection method uses photometry to measure the change 
in received flux from a star when a planet crosses the stellar disc. 
Since then, ground-based transit detection surveys such as WASP 
(Pollacco et al. 2006) and HAT (Bakos et al. 2004) have been set 
up around the world. Once a survey has discovered a transiting ex- 
trasolar planet (TEP), follow-up spectroscopic and photometric ob- 
servations are required to properly characterise the system. 

To achieve high-precision photometry requires not only a high 
signal-to-noise ratio (S/N) but also nullification of many system- 
atic effects inherent in ground-based photometry. Southworth et al. 
(2009a) investigated the use of defocussed telescopes to obtain 
high-precision photometry. Telescope defocussing causes the light 
from point sources to be distributed over thousands of CCD pix- 
els. This allows the use of longer exposure times, which means that 
the CCD is read out less often. This reduced dead time means that 
more photons can be collected over a given time interval, leading 
to lower Poisson and scintillation noise. Flat-fielding noise also av- 
erages down by an order of magnitude due to the large number 
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of pixels in the software aperture, and changes in seeing become 
inconsequential. Many researchers have used the technique of de- 
focussed photometry to obtain precise measurements of the param- 
eters of TEP systems (e.g. Gillon et al. 2007, 2008; Demory et al. 
2007; Winn et al. 2007a,b). We note that the defocussing technique 
is only suitable for bright isolated stars. If a star is too faint then the 
increased background and read-out noise can be detrimental. If the 
target and comparison stars are in a crowded field then defocussing 
will cause their point spread functions (PSF) to become blended 
with those of nearby stars. 

The discovery of the TEP system WASP-50 was presented by 
Gillon et al. (201 1), who found it to comprise a TEP with a mass of 
1.47 ± 0.09 Mj U p and radius of 1.15 ± 0.05 R Jup , orbiting a cool 
star with mass and radius 0.89 ± 0.08 M and 0.84 ± 0.03 R . 
We observed WASP-50 with the aim of improving its measured 
physical properties, using the telescope-defocussing approach. We 
used the New Technology Telescope (NTT) operated by ESO at La 
Silla, Chile. This telescope is alt-az mounted and has a thin primary 
mirror whose shape is maintained using active optics. Whilst these 
qualities are not outwardly well-suited to the telescope-defocussing 
approach, we have previously found the NTT to work well for 
this type of observations (Hellier et al. 2011; Tregloan-Reed et al. 
2013). We observed two complete transits of WASP-50 using the 
NTT, achieving extremely low photometric scatters of 258 and 211 
parts per million (ppm), respectively, versus a fitted model. To our 
knowledge the latter is the lowest scatter ever achieved in ground- 
based photometry per point for a point source. 

Some of the highest photometric precisions previously ac- 
complished for a TEP system are 479 ppm for CoRot-1 using the 
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Figure 1. Surface plot of the PSF of WASP-50 in an image taken at random 
from the observing sequence on the night of 201 1 October 24th. The x and y 
axes are in pixels. The lowest and highest counts are 684 and 24 726 ADUs, 
respectively, and the z axis is on a linear scale. 



8.2 m VLT (Pont et al. 2010), 478 ppm for WASP-4 using the 6.5 m 
Magellan Baade telescope (Winn et al. 2009), 470 ppm for WASP- 
10 using the university of Hawaii 2.2 m telescope (Johnson et al. 
2009), 387 ppm for WASP-2 using a 1.5 m telescope (Southworth 
et al. 2010), and 316 ppm for TrES-2 using the 10.4 m Gran Tele- 
scopio Canarias (Colon et al. 2010). The highest photometric pre- 
cision we are aware of from a ground-based telescope was previ- 
ously 258 ppm in time-series observations of stars in the open clus- 
ter M67 (Gilliland et al. 1993). 

An alternative metric which is well-suited to direct compari- 
son is signal to noise per unit time. We have calculated the scatter 
in ppm per minute of observing time for the datasets listed above. 
By this metric the better of our two datasets is almost exactly equal 
to the best one presented by Gilliland et al. (1993), and both of our 
light curves are better than any previously published ground-based 
photometric observations of a TEP system. 



2 OBSERVATIONS AND DATA REDUCTION 

Two transits of WASP-50 were observed on the nights of 
2011/11/20 and 2011/11/24 using the NTT with the EFOSC2 in- 
strument operated in imaging mode. In this setup the CCD covers 
a field of view of (4.1') 2 with a plate scale of 0.12" px _1 . The im- 
ages were windowed down to 1100 x 1600 pixels and no binning 
was used, resulting in a dead time between consecutive images of 
50 s. The observations were taken through a Gunn r filter (ESO 
filter #784). The moon was below the horizon for half of the first 
transit and all of the second transit. The telescope was initially fo- 
cussed and the shape of its primary mirror was adjusted to obtain 
the best image possible. We then applied a defocus to the telescope 
and performed the full observing sequence without adjusting the 
telescope focus or mirror shape. 

The pointing of the telescope was adjusted to allow a good 
comparison star to be observed simultaneously with WASP-50. The 
comparison star used was 2MASS J02544939- 105 1548, which 
is of a similar apparent magnitude and colour to WASP-50. The 
2MASS J — K 8 colour indices of the two objects are 0.432 for 
WASP-50 and 0.357 for the comparison star (Skrutskie et al. 2006). 
We were able to keep the telescope autoguided through all obser- 
vations. An observing log is given in Table 1. 

Fig. 1 shows the shape of the point spread function (PSF) of 



WASP-50 in an image taken at random from the observing se- 
quence on the night of 2011/10/24. Fig. 1 provides an interest- 
ing comparison with similar observations of WASP-4 and WASP- 
5 from a 1.5 m telescope of more traditional design (see fig. 1 in 
Southworth et al. 2009a and fig. 1 in Southworth et al. 2009b). The 
PSF for the current observations shows a much more rounded an- 
nulus of high counts, which allows a smaller amount of defocussing 
to be used to collect a given number of photons without saturating 
individual pixels. 

We reduced the data in an identical fashion to Southworth 
et al. (2009a). In short, we performed aperture photometry using an 
IDL implementation of DAOPHOT (Stetson 1987), and adjusted the 
aperture sizes to obtain the best results (see Table 1). A differential- 
magnitude light curve was calculated between the target and com- 
parison star. A first-order polynomial was then fitted to the outside- 
transit data and then subtracted to remove a slow trend present in 
the differential magnitudes. The times of the start of the exposures 
were given in JD/UTC in the FITS file headers, and we converted 
these to times of the midpoints of the exposures in BJD/TDB. 

In order to confirm the low scatter of the resulting light curves 
we performed a second data reduction with completely different 
methods. We used the STARLINK/AUTOPHOTOM package (Eaton 
et al. 1999) driven by a custom C-shell script (Southworth et al. 
2004b), and obtained a light curve with an rms scatter of 414 ppm 
for the first night of data. This result agrees with our light curve 
from DAOPHOT, once the discretization of the datapoints (AU- 
TOPHOTOM quotes instrumental magnitudes to only three decimal 
places) is taken into account. 



3 DATA ANALYSIS 

We fitted the WASP-50 data in a similar manner to Tregloan-Reed 
et al. (2013). PRISM 1 (Planetary Retrospective Intergrated Starspot 
Model) was used to model the transit. PRISM uses a pixellation ap- 
proach to represent the star and planet on a two-dimensional array 
in Cartesian coordinates. This makes it possible to model the tran- 
sit, limb darkening and starspots on the stellar disc simultaneously. 
Limb darkening was implemented using the standard quadratic law. 
PRISM used the six parameters given in Table 2 to model the sys- 
tem, where the fractional stellar and planetary radii are defined as 
the absolute radii scaled by the semimajor axis (r SiP = i? SjP /a). 

GEMC 1 (Genetic Evolution Markov Chain) was used to fit the 
model to the data. GEMC begins by randomly generating param- 
eters for the starting points of N chains, within the user-defined 
parameter space, and then simultaneously evolves the chains for X 
generations. At each generation the chains are evaluated for their 
fitness 2 . The parameters of the fittest member undergo a maximum 
±1% perturbation and its fitness is then re-evaluated. If the fitness 
has improved it is accepted but if the fitness has deteriorated it may 
or may not be accepted based on a Gaussian probability distribu- 
tion: 

P = exp ^(xl-D-xU)^ (1) 

where (n — 1) is the previous generational chain and n is the 
current generational chain being evaluated. The next step is to then 
evolve the other chains. This is accomplished in a similar way as 

1 Available from http://www.astro.keele.ac.uk/~jtr 

2 A solution's fitness was found by calculating the 1 /\ 2 value. 



© 0000 RAS, MNRAS 000, 000-000 



Transits in the WASP-50 planetary system 3 



Table 1. Log of the observations presented in this work. N obs is the number of observations. 'Moon ilium.' and 'Moon dist.' are the fractional illumination of 
the Moon, and its distance from WASP-50 in degrees, at the midpoint of the transit. 



Date 


Start time 


End time 


N obs 


Exposure 


Filter 


Airmass 


Moon 


Moon 


Aperture 


Scatter 




(UT) 


(UT) 




time (s) 






ilium. 


dist. 


sizes (px) 


(ppm) 


201 1/1 1/20 


00:59 


06:02 


111 


120-150 


Gunn r 


2.62 -+ 1.48 


0.384 


91.2 


75, 105, 120 


258 


2011/11/24 


01:08 


06:27 


124 


150 


Gunn r 


2.10 1.53 


0.045 


137.4 


75, 100, 125 


211 



a genetic algorithm, in that the chain parameters are modified by 
incorporating the parameters of another chain. But unlike a genetic 
algorithm where a member is picked by a weighted random num- 
ber and then the digits of each parameter are crossed over with the 
digits from a different member, GEMC directly perturbs the param- 
eters of each chain in a vector towards the best-fitting chain. The 
size of this perturbation is between zero and twice the distance to 
the best-fitting chain, allowing the chain to not only move towards 
but also to overshoot the position of the best-fitting chain. This con- 
tinues until all the chains have converged around the optimal solu- 
tion. Once the burn-in phase has been completed the chains cease 
communication and begin independent Markov chain runs. 

Because GEMC is a hybrid between the Markov Chain Monte 
Carlo approach and a genetic algorithm, the burn-in phase is rel- 
atively short, allowing us to use a large parameter search space. 
The boundaries of the search space for each parameter are given 
in Table 2, which also contains the individual results for the two 
light curves. Table 2 also gives the final photometric parameters 
for the WASP-50 system, which are weighted means of the results 
from the two individual fits. All errorbars denote 1-ct uncertainties. 
Fig. 2 and 3 compares the light curves to the best-fitting models, 
including the residuals. The two datasets were modelled individu- 
ally, and the agreement between the best-fit parameters is excep- 
tionally good. The best-fit limb darkening coefficients are also in 
good agreement with the theoretically predicted values for WASP- 
50 A of Ml = 0.407 and u2 = 0.281 (Claret 2004). 

Our data were taken with 120 s and 150 s exposures, so we 
have checked whether these relatively long exposure times affect 
the derived parameters. For this we modelled the data using the 
JKTEBOP code (Southworth et al. 2004a), finding results in good 
agreement with those from PRISM. We then used JKTEBOP's op- 
tion to numerically integrate the model over the duration of each 
exposure whilst finding the best fit (Southworth 2011). The final 
parameters for each light curve altered by only 0.1 to 0.25-a, al- 
lowing us to conclude that smearing of the transit shape due to long 
exposure times does not have a significant effect on our results. 

To check for correlated 'red' noise we used the Monte 
Carlo and residual-permutation algorithms in JKTEBOP (South- 
worth 2008) to assess the uncertainties in the fitted parameters. We 
found a difference between the two methods of only 0.1%, and con- 
clude that correlated noise is not present at a significant level in our 
data. 

We also used JKTEBOP to check whether the removal of the 
slow drift in brightness with a first-order polynomial had any ef- 
fect on our results. We found that including the polynomial coeffi- 
cients as fitted parameters caused changes in the other parameters 
of roughly 0.001 -a. We conclude that the detrending process has 
had no deleterious effect on our results. 

We have collected the available times of mid-transit for 
WASP-50 from the literature (Gillon et al. 201 1; Sada et al. 2012). 
All timings were converted to the BJD/TDB timescale and used to 



Table 3. Times of minimum light of WASP-50 and their residuals versus 
the ephemeris derived in this work. 

References: (1) Gillon et al. (2011); (2) Sada et al. (2012); (3) This work. 



Time of minimum 


Cycle 


Residual 


Reference 


(BJD/TDB - 2400000) 


no. 


(BJD) 




55558.61237 ±0.00020 


0.0 


0.00000 


1 


55849.92131 ±0.00060 


149.0 


-0.00004 


2 


55851.87634 ±0.00028 


150.0 


-0.00010 


2 


55855.78664 ±0.00008 


152.0 


0.00001 


3 


55859.69680 ±0.00012 


154.0 


-0.00001 


3 



Table 4. Physical properties of the WASP-50 system. The equilibrium tem- 
perature, T e ' q , is for an assumed zero albedo and full heat redistribution. 
is the Safronov (1972) number. 



Parameter 


This work 


Gillon etal. (2011) 


M A (M ) 


0.861 ± 0.052 ±0.023 


0.8921°;°?° 


Ra (R©) 


0.855 ± 0.018 ±0.007 


0.843 ±0.031 


log g A (cgs) 


4.509 ± 0.012 ±0.004 


4.537 ± 0.022 


PA (p©) 


1.376 ±0.032 


1 48+ 010 


M b (M Jup ) 


1.437 ± 0.063 ±0.025 


1 AfiR+O- 091 
1.468_ Q 086 


( Rjup) 


1.138 ± 0.024 ±0.010 


1.153 ±0.048 


9h (ms~ 2 ) 


27.50 ±0.64 


27.5 ±1.6 


Pb ( Pjup) 


0.911 ± 0.032 ±0.008 


u.yoo_ 0g2 


^c' q (K) 


1410 ± 26 


1393 ±30 


6 


0.0853 ± 0.0024 ± 0.0007 




a(AU) 


0.02913 ± 0.00059 ± 0.00025 


0.02945 ± 0.00085 


Age (Gyr) 


a i +6.7 +1.5 
B-J- _4.4 -1.3 





obtain an improved orbital ephemeris: 

T = BJD/TDB 2 455 558.61237(20) + 1.9550938(13) x E 

where E represents the cycle count with respect to the reference 
epoch and the bracketed quantities represent the uncertainty in the 
final digit of the preceding number. Fig. 4 and Table 3 show the 
residuals of these times against the ephemeris. We find no evidence 
for transit timing variations in the system. 

3.1 Physical properties of the WASP-50 system 

Once the photometric parameters of WASP-50 had been measured, 
we moved to the determination of its physical characteristics. We 
adopted the approach of Southworth (2009), which uses the param- 
eters measured from the light curves and spectra, plus tabulated 
predictions of several theoretical models. We adopted the values of 
i, r p /r s and r 3 +r p from Table 2, and the stellar properties of effec- 
tive temperature T e g = 5400 ± 100 K, metal abundance = 
-0.12 ± 0.08 and velocity amplitude K s = 256.6 ± 4.4ms" 1 
(Gillon etal. 2011). 
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Table 2. Derived photometric parameters from each lightcurve, plus the interval within which the best fit was searched for using GEMC. 



Parameter 


Symbol 


Search interval 


201 1/1 1/20 


2011/11/24 


Combined photometric parameters 


Radius ratio 


r p /r B 


0.05 to 0.30 


0.13710 ±0.00049 


0.13661 ±0.00036 


0.13678 ±0.00029 


Sum of fractional radii 


r s + r p 


0.10 to 0.50 


0.1552 ±0.0018 


0.1553 ±0.0016 


0.1552 ±0.0012 


Linear LD coefficient 


«1 


0.0 to 1.0 


0.386 ± 0.068 


0.385 ± 0.049 


0.386 ± 0.040 


Quadratic LD coefficient 




0.0 to 1.0 


0.281 ±0.099 


0.279 ± 0.043 


0.280 ± 0.040 


Inclination (degrees) 


i 


70.0 to 90.0 


84.43 ± 0.17 


84.45 ±0.14 


84.44 ±0.11 


Transit epoch (BJD/UTC) 


To 


±0.5 in phase 


2455855.78172 ± 0.000076 


2455859.691755 ± 0.000118 






Figure 2. Transit light curve and the best-fitting model for 201 1/1 1/20. The residuals are displayed at the base of the figure. 





Figure 3. Transit light curve and the best-fitting model for 201 1/1 1/24. The residuals are displayed at the base of the figure. The gap in the data between phases 
0.03 and 0.04 was caused by a technical difficulty with the shutter on EFOSC2. This was corrected by a shutdown and restart of the instrument. The telescope 
pointing was unaffected. 
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Figure 4. Residuals of the available times of mid-transit versus the orbital ephemeris found in this work. The two timings from this work are the last two 
points after cycle number 150. 



An initial value of the velocity amplitude of the planet, K p , 
was used to calculate the physical properties of the system using 
standard formulae and the physical constants listed by Southworth 
(2011). The mass and [^1 of the star were then used as interpo- 
lates within tabulated predictions from stellar theoretical models, in 
order to find the expected T e g and radius. K p was then iteratively 
refined to find the best agreement between the observed and pre- 
dicted T c ff, and the light-curve-derived r B and predicted — . This 
was performed for ages ranging from zero age to the terminal- age 
main sequence, in steps of 0.01 Gyr. The overall best fit was iden- 
tified, yielding estimates of the physical properties of the system 
and the evolutionary age of the star. This procedure was performed 
separately using five different sets of stellar theoretical models (see 
Southworth 2010), and the spread of values for each output param- 
eter was used to determine a systematic error. Statistical errors were 
propagated by perturbing each input parameter in turn to quantify 
the effect on each output parameter. 

The final results of this process have good internal agreement 
(between the five sets of theoretical models) and are also consis- 
tent with those found by Gillon et al. (2011). The final physical 
properties are given in Table 4 and include separate statistical and 
systematic errorbars for those parameters with a dependence on the 
theoretical models. The final statistical errorbar for each parameter 
is the largest of the individual ones from the solutions using each 
of the five different stellar models. The systematic errorbar is the 
largest difference between the mean and the individual values of 
the parameter from the five solutions. 



4 SUMMARY AND DISCUSSION 

In the pursuit of obtaining accurate properties for transiting extra- 
solar planetary systems we have achieved photometric precisions of 
258 and 211 ppm in observations of WASP-50, which to our knowl- 
edge is a record for ground-based photometry of a point source. Our 
approach was to heavily defocus the 3.6 m NTT and to use expo- 
sure times of 120-150 s. We also benefitted from the presence of a 
good comparison star, at a distance of 2.25 arcmin from WASP-50 
and with similar colours and r-band apparent magnitude. The sky 
brightness was also low, as the moon was below the horizon for 
most of our observations. 

We reduced our data using two independent pipelines, finding 
agreement between them. We modelled the light curves using the 
PRISM and JKTEBOP codes, again finding good agreement. From 
these results and published spectroscopic measurements we have 
deduced the physical properties of the WASP-50 system to high 
precision. The properties of the planet WASP-50 b are now known 
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Figure 5. Simulated light curve of a 2 R® planet orbiting WASP-50 A. 



to within 5% (mass), 2% (radius), 4% (density) and 2% (surface 
gravity). This compares to 6%, 4%, 9% and 6%, respectively, in 
the discovery paper (Gillon et al. 2011). We also obtained a re- 
fined orbital ephemeris. Further improvements in precision could 
be made in the shorter term by obtaining additional radial velocity 
measurements, and in the longer term by using sets of stellar mod- 
els which show a better interagreement on properties of the host 
star WASP-50 A. 

In our study of WASP-19 (Tregloan-Reed et al. 2013) we 
found a modest discrepancy between limb darkening coefficients 
measured from three datasets taken with the same telescope. This 
was attributed to the fact that WASP-19 A is an active star with sig- 
nificant starspot activity, which alters the limb-darkening behaviour 
of the star (Ballerini et al. 2012). Whilst WASP-50 does show mod- 
est chromospheric activity, as judged from emission in the Call H 
and K lines (Gillon et al. 2011), starspot anomalies have not been 
observed in any of the five transit light curves of this system. The 
limb darkening coefficents found from our two datasets are in ex- 
cellent agreement (0.02-a), supporting the suggestion that starspots 
affect stellar limb darkening. 

Finally, we have checked if observations of this high precision 
could be used to characterise transiting super-Earths. We injected 
a synthetic transit of a 2R® planet in front of WASP-50 A into 
the residuals of our best fits from both nights, and binned the data 
together. This simulated light curve (Fig. 5) shows a clear transit 
signature, suggesting that ground-based defocussed photometry of 
transiting super-Earths is a viable possibility. 
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